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We show that late Q-ball decay in the MSSM with gauge-mediated SUSY breaking can provide a natural 
source of non-thermal NLSPs which subsequently decay to gravitino dark matter without violating nucleosyn- 
thesis constraints. To show this, we perform a global analysis of Q-ball formation and decay in Affleck-Dine 
baryogenesis for a d = 6 (u c d c d c ) 2 flat direction of the gauge-mediated MSSM. A general phenomenological 
potential for the flat-direction is studied and the Q-ball decay properties are obtained as a function of its param- 
eters. The corresponding gravitino mass necessary to account for dark matter is then determined for the case 
of stau NLSPs. The decay temperature depends on the charge of the Q-balls, which is determined by the frag- 
mentation of the AD condensate. Different fragmentation scenarios are considered, and the final non-thermal 
NLSP density from Q-ball decay and NLSP annihilation is determined. Particular care is taken to establish 
^ 1 that NLSPs from Q-ball decay become homogeneous and non-relativistic prior to annihilation. The gravitino 

, mass necessary for dark matter is naturally consistent with the theoretical gravitino mass in the gauge-mediation 

model. 
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I. INTRODUCTION 



A fundamental question in cosmology is the nature of dark matter (DM) and, in particular, the identity of the dark matter 
particle in the case where it has a particle physics origin. In the context of the MSSM and its extensions, the dark matter 
candidate is fundamentally connected to the nature of SUSY breaking. In gravity-mediated SUSY breaking the neutralino is a 
natural lightest SUSY particle (LSP) candidate. In gauge-mediated SUSY breaking (GMSB), the gravitino becomes a natural 
LSP candidate. 

There are two ways a gravitino DM density might be produced in GMSB models. For sufficiently large reheating temperature 
(N| \ (Tr ~ 10 7 GeV for the case 1113/2 ~ 0.1 GeV), gravitinos can be produced via thermal scattering [1]. Alternatively, gravitinos 
0^ ■ could be produced via the late decay of next-to-lightest supersymmetric particles (NLSPs). The latter process is, however, 
strongly constrained by big-bang nucleosynthesis (BBN) [2]. Decay of conventional thermal relic NLSPs to a DM gravitino 
density is already completely excluded by BBN. Therefore if gravitino dark matter comes from NLSP decay, a source of non- 
thermal NLSPs is necessary. Even in this case, the mass of the gravitino must be below 1 GeV in order to satisfy BBN 
constraints on NLSP decay. The weakest BBN upper bounds are obtained for the case of stau and sneutrino NLSPs, which can 

■ allow gravitino masses up to 0.1-1 GeV 121 • 

A second fundamental question is the origin of the baryon asymmetry. In the MSSM a simple and natural mechanism for 

■ generating the baryon asymmetry is Affleck-Dine baryogenesis [3]. Remarkably, this can also provide a source of non-thermal 
NLSPs. This is because the AD condensate can fragment to Q-balls which subsequently decay to NLSPs at a low temperature 
fill. For the case of d = 6 flat directions, the Q-balls will typically decay to a density of NLSPs at a temperature below 1 GeV, 

' which is well below the NLSP freeze-out temperature 1 . 

Thus, without adding any fields beyond those in the MSSM, it is possible to simultaneously produce both the baryon asym- 
metry and gravitino dark matter via AD baryogenesis. In fact, only Q-ball decay, through the production of a non-thermal 
NLSP density, can produce gravitino dark matter in the case of AD baryogenesis along a d = 6 flat direction 2 . This is because 
the reheating temperature in d = 6 AD baryogenesis must be below about 100 GeV in order to produce the observed baryon 
asymmetry [4], which is too low for thermal generation of gravitino dark matter. 

Before summarizing the analysis to be carried out here, we first state our main result. We will show that it is possible to produce 
gravitino dark matter via decay of the non-thermal NLSPs coming from Q-ball decay. The required dark matter gravitino mass 
can be light enough to be consistent with BBN constraints on NLSP decay to gravitinos. Moreover, the gravitino mass required 
to explain dark matter is naturally consistent with the mass predicted by the GMSB model. Therefore d = 6 AD baryogenesis in 
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' Alternative Q-ball decay scenarios, in which the Q-balls are assumed to decay directly to gravitino dark matter rather than to NLSPs, are presented in 
2 It may also be possible to produce gravitino dark matter via inflaton decay (al or decay of the fields associated with SUSY breaking | 9], but this is highly 
dependent upon the model of inflation and SUSY-breaking. 
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GMSB can produce both gravitino dark matter and the baryon asymmetry without violating BBN constraints. This model can 
therefore simultaneously solve two fundamental problems in the context of MSSM cosmology with GMSB: the origin of the 
baryon asymmetry and the origin of dark matter. 

For the case of a d = 6 flat direction, we showed in [4] that if the messenger mass scale is sufficiently close to the magnitude 
of the AD field <t> at the onset of baryogenesis, AD condensate fragmentation and the resulting Q-balls will typically occur near 
the region of the flat-direction potential between the logarithmic GMSB "plateau" at large |<t>| relative to the messenger mass 
and the quadratic potential at small |<t>|. In this case the Q-balls will be unstable and can decay before BBN. For this to occur, 
the messenger mass scale must be large enough, which in turn requires that j2 is not much smaller than 1 GeV. In this case the 
Q-balls can have properties which are intermediate between the conventional gauge-mediated and gravity-mediated Q-balls 3 . In 
B5D we studied the properties of GMSB Q-ball solutions in detail, in particular how they interpolate between gravity-mediated 
and gauge-mediated type Q-balls and how their decay temperature depends on the potential parameters. 

Our previous studies established the nature of the initial fragmentation of the AD condensate and the form of Q-ball solutions 
in GMSB models. In this paper we consider whether the combined effect of fragmentation to Q-balls and the subsequent decay 
of these Q-balls to NLSPs can explain gravitino DM while remaining consistent with BBN constraints. We will combine the 
semi-analytical fragmentation model of [4] with the numerical Q-ball solutions of [5]. Using these in a global analysis of AD 
baryogenesis, condensate fragmentation and Q-ball formation, we will be able to establish the Q-ball decay temperature to 
NLSPs and so determine the gravitino density and dark matter gravitino mass from the subsequent annihilation and decay of the 
NLSPs. 

Our analysis is based on a phenomenological potential, first introduced in [4], which models a general flat-direction potential 
in the MSSM with gauge-mediation. This potential is characterized by a number of parameters, in particular the soft scalar mass 
m s , the messenger mass M m (above which SUSY breaking is suppressed), the GMSB A-term a and the mass scale of the non- 
renormalizable operators responsible for AD baryogenesis, M. Using this potential, we will evolve the flat-direction condensate 
through AD baryogenesis, condensate fragmentation, Q-ball formation and, finally, decay to NLSPs, in order to establish the 
Q-ball decay temperature Tj and resulting NLSP density as a function of the potential parameters. 

Using the semi-analytical method introduced in J4J], in which fragmentation is studied by perturbing the negative pressure 
AD condensate obtained by averaging over oscillations of the AD field, we will first establish the energy, baryonic charge and 
diameter of the initial fragments. Due to the suppression of the A-terms expected in GMSB models, the energy to (global) 
charge ratio of the initial fragments is much larger than that of a Q-ball of the same charge. The subsequent evolution of such 
fragments has been studied numerically for the case of gravity-mediated Q-balls lllOll . In that case it was shown that the initial 
fragments lose energy by producing pairs of positive and negative charge Q-balls (we will refer to these as '±Q-ball pairs' in 
the following). No numerical analysis of an equivalent resolution exists as yet for the case of gauge-mediation (although earlier 
numerical studies do exist which support the formation of ±Q-ball pairs in gauge-mediation lUllD . However, we can expect a 
similar dissipation of the energy of the initial fragments via the production of ±Q-ball pairs. In the absence of a full numerical 
simulation, we will consider different scenarios for the evolution of the initial fragments to a distribution of positive and negative 
charged Q-balls. To obtain the corresponding Q-ball properties, in particular their decay temperature, we must then match the 
energy of the initial fragments to the energy of the corresponding final state GMSB Q-balls, following the methods of ||5|]. 

Given the Q-ball decay temperature and the nature of the NLSP, the final NLSP density from Q-ball decay can then be 
obtained. In order to do this, it is important to consider in detail how Q-balls decay and how the initially relativistic NLSPs from 
Q-ball decay lose energy and disperse. Q-balls are spatially separated when they decay, therefore if the NLSPs were to rapidly 
stop via scattering with thermal background particles and annihilate while inhomogeneous, a much lower NLSP density would 
remain than in the case where the NLSPs can homogenize prior to annihilation. We will show explicitly that the NLSPs from 
Q-ball decay generally homogenize via relativistic free-streaming well before they annihilate. Therefore the final NLSP density 
is determined by non-relativistic annihilation of a homogeneous NLSP density produced by Q-ball decay, as is conventionally 
assumed. 

The NLSPs subsequently decay to gravitinos. If the gravitinos account for the observed dark matter density, they should have 
a mass less than 1 GeV in order that NLSP decay is compatible with BBN [2]. In the case where only positively charged Q-balls 
form in AD baryogenesis, this is a serious obstacle to producing gravitino dark matter via Q-ball decay, since B-conservation 
combined with R-parity conservation implies that the dark matter particle must satisfy iuom ~ 2 GeV HlUdUO"!!] 4 . However, 
in the case where almost neutral initial fragments break up into positive and negative charged Q-balls, a larger NLSP density 
will be produced from decay of the ±Q-ball pairs. This is a key feature of the GMSB Q-ball decay model for gravitino dark 
matter. For the case of stau NLSPs, we will show that the subsequent annihilation of the NLSP density can naturally produce 



3 By conventional gauge-mediated Q-balls, we mean Q-balls in which |<J>| is large enough that the solutions are completely dominated by the logarithmic 
potential. 

4 This is also an obstacle to dark matter in the gravity-mediated MSSM. A possible solution in that context was proposed in Qj] , in which the LSP is an axino 
of mass f» 2 GeV. 
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the observed dark matter density with gravitinos in the mass range 0. 1-1 GeV. 

Finally, by considering the simplest model for GMSB, based on a vector pair of SU(5) messenger supermultiplets, we will 
show that the gravitino mass required for dark matter is naturally compatible with the corresponding theoretical gravitino mass 
from GMSB. 

Our paper is organized as follows. In Section 2 we review the flat-direction potential and the semi-analytical method of [4] 
for determining the nature of the initial AD condensate fragments. We also introduce a simple model for the evolution of these 
fragments into pairs of positive and negative charged Q-balls. In Section 3 we discuss GMSB Q-ball solutions and compute 
the Q-ball decay temperature to NLSPs. We also establish the conditions under which the NLSPs become non-relativistic and 
homogenize prior to their annihilation. In Section 4 we discuss annihilation of the NLSPs from Q-ball decay and establish the 
gravitino mass necessary to account for dark matter as a function of the Q-ball decay temperature T £ /. We also discuss the BBN 
constraints on NLSP decay. In Section 5 we present the results of our global analysis of AD baryogenesis and gravitino dark 
matter via Q-ball decay. We follow the evolution of flat-direction field from the inflaton-dominated era through the generation 
of the baryon asymmetry, the fragmentation of the condensate to an ensemble of ±Q-ball pairs and their subsequent decay to a 
gravitino density via NLSP annihilations and decays. We determine the Q-ball decay temperature and required gravitino dark 
matter mass as a function of the flat-direction potential parameters. In Section 6 we compare the gravitino mass necessary to 
account for dark matter with the theoretical gravitino mass from GMSB. In Section 7 we present and discuss our conclusions. 



II. FLAT-DIRECTION POTENTIAL, AFFLECK-DINE BARYOGENESIS AND CONDENSATE FRAGMENTATION IN GMSB 

A. GMSB Flat-direction Potential 

In iU we introduced a phenomenological potential which models the expected behaviour of a GMSB flat direction potential 
as a function of the flat direction field <t>. We will focus on a d = 6 (u c d c d c ) 2 flat direction in the following. The potential is 



V(4>)=mX,ln 2 I J + m 2 /2 ^1 +*ln ) ) |4f - cH^ + (AW +h.c.) - 

Here 



4>5 2 



5!M 3 



(1) 



(1)6 

W = oW (2) 

is the effective superpotential for the flat direction superfield <t>. We include the factor 6! in Eq. (f2]) so that the physical strength 
of the interactions is dimensionally of the order of M. The scale M is usually assumed to be of the order of the Planck scale, 
however we will consider it to be a free parameter throughout this paper. The log-squared factor in first term in Eq. (fl} is due to 
GMSB with messenger mass M m lfl5ll . The factor multiplying this takes into account 1-loop radiative corrections due to gauging 
loops once g|<t>| < M m , with K ps —(0.01 — 0.1) (analogous to gravity-mediated SUSY breaking corrections lfl2l IT3I flql - ). The 
Standard Model (SM) couplings g at the renormalization scales of interest (p ~ M m ~ 10 13 GeV) are g s» 0.6 — 1; for simplicity 
we will set the factor g to 1 in V(<P). The second term is due to gravity-mediated SUSY breaking including the 1-loop correction 
term K; for simplicity we will set K = K. For the A-term we consider 

a a ° ms ri\ 
A = m V2 + - ^T72- (3) 

1 + M 

The first term in Eq. © represents the A-term due to gravity-mediated SUSY breaking. The second term models the A-term in 
gauge-mediated SUSY breaking at |4>| < M m , which is generated at 1-loop from the gaugino masses and is therefore suppressed 
compared with the A-term in gravity-mediated models. We will consider a u = a/ (4k) « 0.01 to be a typical value, where 
a = g 2 /4K and g is the relevant SM gauge coupling. (We will however consider the effect of varying a .) The suppression factor 
(1 + |4>| 2 /M 2 )-'/2 models the 1 /|4>| suppression of the GMSB A-term at |4>| > M m fiUl . 



B. Affleck-Dine Baryogenesis 

As H decreases, the |<t>| 7^ minimum becomes unstable and the <t> field begins oscillation in its real and imaginary directions. 
The B-violating A-terms induce a phase difference between the oscillations in the real and imaginary directions, such that at late 
times a 3 / 2< t> (where a 3 / 2 takes into account the effect of expansion) describes an ellipse in the complex plane, corresponding 
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to a baryon asymmetry Due to the suppression of the A-terms in GMSB, the condensate is typically highly elliptical (a 
precessing ellipse, to be precise (01). This means that there is much more energy in the GMSB AD condensate than would be 
expected if it were entirely made of baryonic charge (corresponding to a circular condensate, so-called "Q-matter"), therefore 
the condensate is largely neutral with respect to global charge. This fact will play a key role in the Q-ball decay scenario for 
gravitino DM. 

AD baryogenesis occurs during matter domination by inflaton oscillations. The degree to which the baryon density is sub- 
sequently diluted relative to the entropy density depends on the duration of inflaton matter domination, with greater dilution 
corresponding to a lower reheating temperature, Tr. Tr is chosen to reproduce the observed baryon asymmetry. 



C. Initial Fragmentation of the AD Condensate 

Once the baryon asymmetry is established, perturbations of the AD condensate, corresponding to primordial density pertur- 
bations, will grow. The condensate will subsequently fragment when these perturbations are comparable to the mean energy 
in the condensate. In ]4J] this was studied by following the growth of perturbations due to the average negative pressure in the 
condensate, which allows the expansion rate at fragmentation, Hf rag , and the wavelength of the dominant perturbation mode, 
Xfrag, to be determined for a given set of potential parameters. This then determines the energy and baryonic charge of the initial 
AD condensate lumps. 

We first summarize the method of JU]. The key observation is that the average (negative) pressure in an elliptical condensate 
is equal to the pressure in a circular condensate, with the constant |<t>| of the circular condensate set equal to the amplitude 
of oscillation of the elliptical condensate. In this way the rate of growth of perturbations oc(f ) of an elliptical condensate can 
be studied using the same method as for a circular condensate [ 17]. This greatly simplifies the analysis of fragmentation as a 
function of the model parameters, since it allows a semi-analytical approach which enables a wide range of model parameters 
to be easily scanned. The semi-analytical approach also has the advantage of being robust, in that it is independent of numerical 
issues such as the resolution of a lattice simulation. It can therefore complement numerical simulations. 

The growth of perturbations of the amplitude and phase satisfies Sep, 80 e s ^\ where 4> = q>e' e /V2. The rate of growth is 
then a = S, where HO 

2 |k| 2 1 (-7 " |k| 2 6 4 \ 
a 2 (V" + 3i2 2 ) ^ a 2 (y* +3&) 2 J 

Perturbation growth is possible once the physical wavenumber \k/a\ of a given perturbation is smaller than a maximum value, 
%nax/a\, where 

= K^(J°L)x(tf-V"). (5) 



V" + 3Q. 2 

The growth of a perturbation of (comoving) wavenumber k, once it is smaller than the critical wavenumber for growth, is then 
given by 



S(t) = f a{t)dt , 



(6) 



where f» is the time at which the mode begins to grow, corresponding to the time when |k/a| = K m . The subsequent growth of 
the perturbation follows from 



S(t) = / 

JO, 



a a(k,a) 



aH 



da , (7) 



where a* is the scale factor at and 



Ikl 2 



1/2 



a(k,a) re M-tH; ■ (8) 

1 - 2 (V +3£2 2 ) 1 

With H °c a~ 3 / 2 during inflaton-domination, this gives 

S(k,a) = 2a(k,a) ( 1 - I H l . (9) 




5 



At a given value of a, the mode with the maximum growth is found by maximizing S(k,a) with respect to |k| . With a(k, a) °c |k| 
and a*(k) = K~ |k|, this is maximized at 



^ = ^ . (10) 
a 9 

The value of S for this mode is 

S=^aH- 1 . (11) 

Here a is determined by the value of Q. and V for the condensate at a given a. The condition for the condensate to fragment is 
that 

SR 8R c . 

Y = ~~R~f ~ 1 ' (12) 

for the mode with wavenumber Eq. ( TTOb . 8R /R ss 5p/p « 10~ 4 is due to the primordial density perturbation, which also 
perturbs the energy density of the condensate. The value of |k| mul /a when this is first satisfied, \k\f rag /a, will determine the size 
of the fragments. The diameter of the fragments is then given by 

fyrag ~ "7 T • (13) 

(J k l/™i>J 

Their global 1/(1) charge (Q = 3B) is 

Q~^fra g nQ (14) 

and the energy of the initial condensate fragment is 

E~tf rag p, (15) 
where we consider each fragment to come from a cube of side "kfrag- 



D. Evolution of the initial condensate fragments 

The initial fragments are almost neutral, corresponding to a highly elliptical condensate. Therefore, the energy to (global) 
charge ratio, E/Q, of the condensate fragments in GMSB models is large compared with the AD scalar mass. As a result, the 
fragments will lose energy to reach a lower energy state. The scalar field dynamics of this process is highly non-linear and 
can only be fully understood via a numerical simulation of the evolution of the fields. A high-resolution simulation has been 
performed for the case of AD baryogenesis in gravity-mediated SUSY breaking models [10], but as yet there is no equivalent 
simulation for gauge-mediated models 5 . Nevertheless, since the non-linear dynamics is expected to be similar in the gauge- 
mediated case (in particular, in the limit where the AD field at initial fragmentation is close to or less than the messenger mass, 
in which case the potential will be similar to a gravity-mediated potential), we can use the results of ifioll as the basis for a simple 
model for the evolution of the condensate fragments in gauge-mediated models. 

The results of H 1 CW1 can be summarized as follows. They consider an AD condensate forming initial fragments with E /mQ « 
1 /e. (m is the AD scalar mass in 1 1011 .) The case of most relevance is a 3D simulation with e = 0.01 (denoted "3D5" in [ 10]). 
This is a gravity-mediated flat-direction with K = —0.1. In this case the initial condensate fragments from AD baryogenesis 
(called "first generation Q-balls" in fToll ) break up into a distribution of ±Q-ball pairs with E/m\Q\ w 1 . The charge distribution 
of the Q-balls (Figure 25 of llOfl) shows that the mean charge \Q\ of the final Q-balls is approximately 4-5 times larger than 
the mean charge of the initial fragments. The initial fragment energy is Ef mg w 100m\Q\f rag , where |<2|/rag is the charge of the 
initial fragment. Therefore with final charge of the Q-balls given by \Q\ w (4 — 5)|Q|/ m? we find Ef rag w (20 — 25)m|g|. The 
energy of the final Q-balls is Eq ss m\Q\. Therefore the ratio of the initial fragment energy to the final Q-ball energy in the 3D5 
simulation of iflbTl is Ef rag s» (20 — 25)£g. 



5 Earlier GMSB simulations do exist 1 1 1], which support the break-up of the initial fragments into ±Q-ball pairs. 
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Based on this, we will model the decay of the fragments to ±Q-ball pairs by a simple model in which the initial fragments 
break up into n p pairs of ±Q-balls of equal energy. (The complete fragmentation process is illustrated schematically in Figure 
1.) The energy of the Q-balls is then Eg = Ef rag /2n„. The limiting case corresponds to n p — 1, with the fragments breaking 
up into a single pair of ±Q-balls. The simulation of jlOll indicates n p « 10—15 may be more realistic, therefore we will also 
consider larger n p . 

Affleck-Dine Condensate breaks up 
into initial nearly neutral fragments 

/ \ 



/ \ 

Neutral fragments decay 
QQ into n p + Q-ball pairs qQ 

f0 90 




Clusters of Q-balls from 
initial fragments decay to NLSPs 
at T d when separated by distance d s 



dplgo* , ,oi 

° d 

s 



FIG. 1: Schematic evolution of the Affleck-Dine condensate fragments. The nearly neutral initial condensate decay into pairs of positive (light) 
and negative (dark) Q-balls. The resulting clusters of Q-ball later decay to NLSPs at when separated by d s , The individual Q-balls in the 
clusters are separated by d p at Tj. 



III. GMSB Q-BALL SOLUTIONS AND Q-BALL DECAY 



The squarks which form the Q-ball decay primarily to quarks and gauginos, with the gauginos subsequently decaying to 
NLSPs. (If decay to gauginos is kinematically excluded then the Q-balls decay directly to quarks and NLSPs.) The Q-ball decay 
rate to fermions has been calculated in [18]. In the saturated case, which occurs when g<p 3> CO, where <p„ is the magnitude of 
the field inside the Q-ball and CO is the effective mass of the Q-ball squarks, the decay rate of the Q-ball is given by Mlal 

CO 3 /? 2 

r "=48^< (16) 

where R is the radius of the Q-ball. The Q-balls decay once T ( / ss H, with the expansion rate given by H = kjT 2 /M p , where 
kj = (iz 2 g(T) /90) x l 2 and g(T) is the effective number of massless degrees of freedom. Therefore the Q-ball decay temperature 
is 

t(lfR 2 M^ 'I 1 
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In order to estimate the Q-ball decay temperature, we therefore need to obtain the Q-ball parameters co and R and the Q-ball 
charge Q of the ±Q-balls. 

Our method is to find GMSB Q-ball solutions which match the Q-ball energy from the fragmentation model, Ef rag /2n p , for 
a given n p . The GMSB Q-balls are minimum energy solutions for a fixed global charge Q. The potential is given by Eq. (fl}, 
but neglecting the A-term, H 2 mass term and non-renormalizable term, which are negligibly small inside the Q-ball. The Q-ball 
solutions have the form 



<F 



9(r) 

7f 



where the Q-ball profile <p(r) is given by the solution of 



3 2 cp i 2 3ip _ dV 2 

r dr 8(p ' 



dr 2 



(18) 



(19) 



with boundary conditions cp (r) = as r — >• and (p(r) -^Oasr->=o. The total charge and energy of the Q-ball are then 



Q= f 4%r 2 (oy(r) 2 dr ; E = [ \%r 2 
Jo Jo 



l/9cp 
2 I 3r 



-y((p) 



co 2 (p 2 



dr. 



(20) 



For each value co there is a unique Q-ball solution with energy E and charge Q. By matching the energy E from the fragmentation 
model to the corresponding Q-ball solution with the same energy, the values of CO, Q, E/Q and R can be obtained 6 . 

The analysis of [18] is for a thin-walled Q-ball of radius R, with a step-like profile for <p(r). In reality, Q-balls in the MSSM 
are not generally thin-walled. However, in the limit where the field in the Q-ball significantly larger than the messenger mass, 
the Q-ball has a broad profile [5] and a thin-walled approximation is a reasonable approximation to the Q-ball. In the opposite 
limit, where the field is small compared with the messenger mass, the Q-balls will be similar to gravity-mediated type Q-balls, 
which are thick- walled with a Gaussian-like profile. A very recent analysis [19] has considered Q-ball decay with a general 
profile. It is shown that, in the saturated limit, the decay rate is similar to Eq. ( fToT ) for the case of GMSB Q-balls with a large 
field strength compared to the messenger mass and somewhat larger that Eq. (fT6b (typically by a factor ~ 20) for the case of 
gravity-mediated-like Q-balls. In the latter case the decay temperature can be a factor ^4 — 5 larger than Eq. ( fT7b . In fact, we 
will show that the Q-balls with field strength large compared with the messenger mass are the ones primarily responsible for 
gravitino dark matter which is consistent with BBN, therefore Eq. ( fT6] l should be a good approximation to the Q-ball decay rate. 



IV. NLSPS AND GRAVITINO DARK MATTER FROM Q-BALL DECAY 
A. Q-ball decay to NLSPs 

Provided that the mass of the NLSP, m x , is less than E/Q, decay to NLSPs is kinematically possible. The NLSPs will initially 
be produced with relativistic or near-relativistic energy 0(E/Q). The final NLSP density will then depend on the subsequent 
evolution of the initially relativistic NLSPs. 

The conventional picture of Q-ball decay assumes that the NLSPs homogenize and become non-relativistic prior to annihi- 
lation lf20l I2H1 . It is important to check this assumption. The NLSPs from Q-ball decay will lose energy via scattering with 
relativistic particles in the thermal background. They will then slow and become non-relativistic within a distance d^R. At 
Q-ball decay, the Q-balls from different initial fragments will be in clusters separated by a distance d$ (Figure 1). Therefore 
if dm 3> ds, the NLSPs will cross between the Q-ball clusters via relativistic free-streaming before they slow by scattering 
and become non-relativistic. In this case the NLSPs will become homogeneous and non-relativistic before they annihilate. We 
discuss the calculation of d^R and ds in the Appendix. We will show in the next section that d^R 3> d$ is easily satisfied for the 
GMSB Q-balls of interest to us here. We can therefore compute the final NLSP density under the assumption of a homogeneous 
and non-relativistic NLSP density from Q-ball decay. 

We next consider the evolution of the NLSP density from Q-ball decay. We first consider a density of neutral NLSPs. Q-balls 
decay away completely over a time H {T4). From Eq. ( fT6b the Q-balls decay according to 

dQ co 3 /? 2 

-~ = -K d ; K d = . (21) 

dt 4871 V 



6 Since in general the Q-ball profiles are not exactly step-like, as assumed in flUl . we will define R in our solutions to be the radius within which 90% of the 
total energy is found. 
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Therefore the homogeneous non-relativistic NLSP density n will increase at a constant rate as the Q-balls decay and inject 
NLSPs. Eventually the annihilation rate of the NLSP density, T ann =n <av>, where o is the annihilation rate and v the relative 
velocity of the NLSPs, and the average is over the NLSP distribution, will become equal to the expansion rate, at which point 
annihilations become effective and the density no longer increases. The final density of non-thermal NLSPs from Q-ball decay 
is therefore 

n M m{Td) « ~~~~ • ( 22 ) 

This density will be achieved provided that the total NLSP density from Q-ball decay in the absence of annihilations, «ror(7d), 
is larger than n ann {T ( [). This will be shown to be generally true in the cases of interest to us here. 

In the case where the NLSPs carry a gauge or global charge, as in the case of stau NLSPs, the Q-balls will produce equal 
number densities of NLSPs and anti-NLSPs. In this case n ann (Td) in Eq. ( 1221 gives the number density of the NLSPs (which is 
equal to the number density of anti-NLSPs), the total NLSP number density is the sum of NLSP and anti-NLSP number densities 
and < ov > in Eq. d22l) is the total annihilation cross-section of an NLSP with both NLSPs and anti-NLSPs from Q-ball decay. 



B. Dark matter gravitino mass 

Given the NLSP density from Q-ball decay, we can determine the value of the gravitino mass m necessary to account for 
the observed dark matter density. By R-parity conservation, one gravitino is produced per NLSP decay, therefore = n NLSP- 
The gravitino density at present is then 

n3l2{Ty) = ^&n ann {T d ). (23) 
8(Td) T d 

Therefore gravitino mass required to account for a dark matter abundance ^3/2 is W3/2 = 9^3/2/^/2, where p 6 is the critical 
density. 

We next apply this to the case of stau NLSPs, which we will focus on in the following. There are two dominant annihilation 
channels for a stau 1 2211 : (i) XX* — !• yy, 7Z(where x* denotes the anti-stau), and (ii) vt — > XT. The average cross-sections times 
relative velocity for these processes are i22ll 

<CTV>W .«^» (24) 



I6na 2 em m 2 s 

< ov >«« s—-^ , (25) 

c 4 Q / ™2 — 1 m 2 



and 



cos 4 B w [ m~ + nig 

where B is the bino. The stau annihilation rate is then 

r fl „„ = «x < ov >xx +n=t* < ov >xx* • (26) 

Therefore, using n T = n^*, we have 

r fl „ n = n=t<av> ; < ov >=< ov >^ + < ov • (27) 
The number of staus from Q-ball decay is then 

< ov > 

The same number of x* are produced. The total NLSP density is therefore 



(28) 



, T s . _ 2H(T d ) 

nNLSP (Td ) = m + nv . (29) 

< Ov > 
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m 3/2 ~ 0.06 GeV x (f ^ ^ , ( _ , . (31) 



From Eq. ( f24b and Eq. d25l l. we see that < ov > xx <C < Ov > xx * in the limit where nig ^> ra x , while < ov > xx and < ov > xx * 
become comparable at mg ps ra x . In the following will assume < ov > xx =< Ov > xx *, noting that < ov > could be smaller by a 
factor of 2 if the bino is sufficiently heavy. 

Using < ov > = 2 < ov > xx *, we find that the gravitino mass required to account for the observed density of dark matter in 
the case of stau NLSPs is 

_ g(T d ) T d Cl 3/2 p c Ana 2 em M p 
g{Ty) krjj mi 

Therefore 

'g(T d )\ l/2 { T d \(\TfN\ 2 

In this have used g(T d ) = 10.75 for T d ~ 10 MeV, a em = 1/128, £ly 2 h 2 = 0.1 126 S and p c = 8.1/z 2 x 10~ 47 GeV 4 . 

C. BBN Constraints 

The BBN constraints on NLSP decay to gravitinos have been analyzed in detail in J3l- The constraints follow from the decay 
of NLSPs at T < Tbbn ~ 1 MeV, which can modify light element abundances formed during BBN. The results of 10] can be 
summarized as follows. Non-thermal stau or sneutrino NLSPs can produce gravitino dark matter and remain consistent with 
BBN if the NLSP mass is greater than around 300 GeV and my 2 < 1 GeV (Figs. 14 and 16 of 10]), while a bino NLSP of mass 
greater than about 300 GeV is consistent with BBN if 

my 2 ~2x 10~ 2 GeV (Figs. 9 and 10 of H). 

There is a large hadronic 

component in the decay products of bino decay, which leads to stronger constraints than the case of a stau or sneutrino NLSP, 
which have primarily radiative decays with only a small hadronic component. (The stau can also catalyze formation of 6 Li via 
formation of 4 He-T bound states.) The analysis of |0] considers only two NLSP masses, 100 GeV and 300 GeV. However, the 
upper bound is generally expected to become weaker for larger NLSP mass, due to the earlier time of decay. 

The upper bound on the gravitino mass can be understood as a lower bound on the time of NLSP decay. The NLSP decay 
rate is proportional to rn 5 NSp /my 2 . Therefore, for a given NLSP mass, the NLSPs decay earlier as the gravitino mass becomes 
smaller. If the NLSP has primarily radiative decays, as in the case of stau and sneutrino NLSPs, then as long as the decay occurs 
sufficiently before the formation of light element abundances at T ps 0.08 MeV, the electromagnetic cascade from NLSP decay 
will not dissociate the light elements. 

Therefore, in order to evade BBN constraints on stau decay to gravitinos, the dark matter gravitino mass must be less than 1 
GeV. From Eq. (ED we see that this requires that the Q-balls decay at a temperature of the order of 10 MeV. 



V. GLOBAL ANALYSIS OF AFFLECK-DINE CONDENSATE FRAGMENTATION, Q-BALL DECAY AND GRAVITINO DARK 

MATTER 

In this section we will present the results of our global analysis of AD baryogenesis, Q-ball formation and decay to gravitino 
dark matter, using the methods discussed in the previous sections. We will focus on the case of stau NLSPs. 

There are a number of conditions that must be satisfied in order to have successful production of gravitino dark matter via 
Q-ball decay to stau NLSPs: (i) nucleosynthesis requires that T d > 1 MeV, in order that the baryon asymmetry is released from 
the Q-balls before the onset of nucleosynthesis at Tbbn ~ 1 MeV; (ii) the gravitino mass should be less than 1 GeV, in order that 
stau decay does not alter element abundances from BBN; (iii) the Q-balls should be able to decay directly to stau NLSPs. This 
requires that E/Q > m%. The present LHC CMS lower bound on the stau mass is ra x > 223 GeV [24]. We will conservatively 
set m x = 250 GeV throughout; larger values of m x will give lower values of tnyi. 

In Table 1 we present the Q-ball decay temperature and dark matter gravitino mass as a function of m s for the case n p = 10, 
a = 0.01, M = M p and M m = 10 13 GeV. (We set K = —0.01 throughout; our main results are not sensitive to the value 
of K, which only becomes important for gravity-mediated type Q-balls.) We consider large values for the AD squark mass, 
m s > 1 TeV, as suggested by the non-observation of squarks at the LHC. We find that dark matter gravitino production satisfies 
all BBN and experimental constraints in all the examples shown in Table 1 . This is our main result: gravitino dark matter can 
successfully originate from Q-ball decay to NLSPs in GMSB. 

The gravitino mass for the cases in Table 1 is in the range 0.1-1 GeV, while the Q-ball decay temperature is in the range 
4.4-7.7 MeV. (In the case where the bino mass is much larger than the stau mass, my 2 can be smaller by a factor of about 2 
for a given T d .) The value of E/Q is in the range 400-706 GeV. The Q-balls generally have cp(0)/M m > 10, therefore they are 
gauge-mediated type Q-balls, dominated by the logarithmic plateau of the potential. 
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In Table 2 we show the corresponding d = 6 GMSB AD baryogenesis and condensate fragmentation parameters. The reheating 
temperature is chosen to reproduce the observed baryon asymmetry and is in the range 44-87 GeV. 

In Table 3 we give the values of d^R and ds for the Q-balls in Table 1. We find that dux ^> ds, with dun ~ 10 18 GeV -1 and 
ds ~ 10 13 GeV -1 . Therefore the stau NLSPs from the decay of these Q-balls will homogenize by free-streaming and become 
non-relativistic well before annihilating. This is generally true of all the cases we consider here. 

In Tables 4 to 7 we show how these results depend on the other parameters of the model. We show the effect of varying M in 
Table 4, the effect of varying M m in Table 5, the effect of varying a a in Table 6 and the effect of varying n p in Table 7. All tables 
use the first entry of Table 1 as their starting point. 

From Table 4 we find that larger values of M are preferred. Decreasing M from M p causes mi, n to rapidly increase above 1 
GeV. Thus the value of M which is often considered the most natural value, M ~ M p , is also preferred by gravitino dark matter 
from Q-ball decay. 

From Table 5 we find that decreasing the messenger mass M m causes both T ( ( and the gravitino mass to decrease. There is 
a small window fromM,,, = 8 x 10 12 GeV to 1 x 10 13 GeV for which Tj is greater than 1 MeV while m 3 / 2 is less than 1 GeV. 
Larger M m causes the gravitino mass to violate the BBN bound while smaller M m causes Tj to violate the BBN bound. The 
gravitino mass can be decreased significantly by considering either a larger bino mass or a stau mass larger than 250 GeV. For 
example, increasing the stau mass to 500 GeV would decrease the gravitino mass for the case M m = 2 X 10 13 GeV in Table 5 to 
0.63 GeV, so increasing the upper bound on the allowed range of M m values. 

From Tables 4 and 5 and Figure 2 we find that there is a strong correlation between the value of <p(0) /M m , where (p(0) is the 
value of the field at the centre of the Q-ball, and the dark matter gravitino mass. The Q-ball decay temperature and gravitino mass 
rapidly increase as (p(0) /M m approaches 1. This is due to the transition of the Q-balls from gauge-mediated to gravity-mediated 
type. This is confirmed by the fact that CO — > m s in the limit where <p(0) /M m <C 1, as expected for gravity-mediated type Q-balls. 
We also find the <p(0)/M„, and the gravitino mass become very sensitive to the model parameters once <p(0)/M„, is significantly 
larger than 1. This is due to the insensitivity of the potential to <p at large tp/M,„, where the potential has only a log-squared 
dependence on <p. As a result, small changes in the potential energy can produce large changes in the value of the field. 

From Table 6 we find that the Q-ball decay temperature and T c / is essentially independent of the GMSB A-term a . At first 
sight this seems surprising. However, in GMSB models the baryon asymmetry contributes only a small part of the charge of the 
± Q-balls coming from the almost neutral initial fragments. Therefore the a a term has very little effect on the Q-balls, provided 
that the AD condensate is strongly elliptical. In this case the only effect of the a term is to alter the net baryon asymmetry in 
the AD condensate and so the reheating temperature, however this plays no role in Q-ball decay. 

From Table 7 we find that decreasing n p , corresponding to production of larger charge ±Q-ball pairs from the initial conden- 
sate fragments, decreases both and m^/i- Thus smaller n p is preferred by the gravitino BBN bound, provided that the Q-ball 
decay temperature is greater than 1 MeV. 

In obtaining these results we have assumed that a sufficient number density of NLSPs is produced from Q-ball decay to be 
able to reach the annihilation density riNisp{T c i) in Eq. d29l , This requires that rtjor > nNLSp{Td), where hjot is the total NLSP 
number density from complete Q-ball decay in the absence of annihilations, utot is given by 

"7OT ~ —jr- ■ (32) 
d s 



This follows because the Q-ball clusters are separated by ds at decay and contain 2n p Q-balls of charge ±Q, with one NLSP 
being produced per unit Q-ball charge. 

For the first entry in Tables 1-3, we have n p = 10, ds = 1.4 x 10 13 GeV -1 and Q = 5.6 x 10 24 , which gives hjot = 4.1 x 
10~ 14 GeV 3 . With T d = 4.4 MeV, k Td ~ 1 and m~ x = 250 GeV, Eq. ([29]) gives n NLSP (T d ) = 6.5 x 10~ 16 GeV 3 . Therefore n TOT > 
nNLSp{Td) and so the NLSPs will annihilate. This also makes clear the importance of annihilations in limiting the NLSP density 
from Q-ball decay. In the absence of annihilations, the NLSP density would be ~ 10 2 times larger, requiring a very much smaller 
gravitino mass. This would not be compatible with the large messenger mass necessary to keep E /Q large enough to allow decay 
to NLSPs. 



m s ( TeV) 


M( GeV) 


M m ( GeV) 


(p(0) /M m 


C0 2 /ra 2 


E( GeV) 


Q 


E/Q{ GeV) 


T d { GeV) 


m 3/2 (GeV) 


1.0 


2.43£18 


1£13 


15.2 


0.108 


2.23£27 


5.59£24 


399 


0.0044 


0.40 


1.5 


2.43£18 


1£13 


17.4 


0.093 


1.97£27 


6.45£25 


557 


0.0067 


0.55 


2.0 


2.43£18 


1£13 


18.9 


0.083 


2.07£27 


2.79£24 


706 


0.0077 


0.67 



TABLE I: Q-ball properties and dark matter gravitino mass for n p = 10 and a — 0.01 as a function of varying m s . 
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m s (TeV) 


M( GeV) 


M m (GeV) 


T R { GeV) 


H fmg (GcV) 


^^(GeV- 1 ) 


|*| /raK /M m 


Efragi GeV) 


Qfrag 


1.0 


2.43£18 


1E13 


44.0 


4.27 


0.050 


2.60 


4.46£28 


1.17E23 


1.5 


2.43£18 


1E13 


66.1 


6.07 


0.035 


2.85 


4.17£28 


5.43£22 


2.0 


2.43£18 


1E13 


86.9 


7.69 


0.027 


3.04 


3.97£28 


3.13£22 



TABLE II: AD baryogenesis and initial condensate fragmentation parameters for n p = 10 and a„ = 0.01 as a function of varying m s . 



m s ( TeV) 


M( GeV) 


M m (GeV) 


d s ( GeV" 1 ) 


^(GeV- 1 ) 


1.0 


2.43£18 


1£13 


1.46£13 


2.98£18 


1.5 


2.43£18 


1E13 


8.53£12 


9.62£17 


2.0 


2.43£18 


1£13 


5.77£12 


4.31E17 



TABLE III: Q-ball decay parameters for n„ = 10 and a„ = 0.01 as a function of varying m s . 



VI. COMPARISON OF THE GRAVITINO DARK MATTER MASS WITH THE THEORETICAL GMSB GRAVITINO MASS 

In the previous section we showed that it is possible for Q-ball decay in GMSB models to produce non-thermal NLSPs which 
have the right density to account gravitino dark matter of mass less than 1 GeV. However, a given GMSB model will also predict 
the gravitino mass. We therefore must compare the dark matter gravitino mass with the theoretical gravitino mass. 

We will consider a standard GMSB framework with n$ vector pairs of messenger chiral supermultiplets f+f, each transform- 
ing as a 5 + 5 of SU(5) j25ll . The superpotential coupling of the messenger supermultiplets to the SUSY breaking supermultiplet 
is 

W = KSff . (33) 

The S supermultiplet has expectation value S —< S > +Q 2 < F$ >. In this case the masses of the / and / scalars and fermions 
are 




5e17 1e18 1.5e18 2e18 2.5e18 



M (GeV) 

FIG. 2: Plot of m 3 / 2 (dashed line, GeV units) and 10(p(0)/M m (solid line) versus M for the case M m = 10° GeV. As (p(0)/M„, approaches 1, 
corresponding to the transition from gauge- to gravity-mediated type Q-balls, the dark matter gravitino mass rapidly increases. 
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m s ( TeV) 


M( GeV) 


M,„( GeV) 


(p(0) /M,„ 


CO 2 /m] 


E{ GeV) 


Q 


E/Q( GeV) 


T d { GeV) 


m 3/2 (GeV) 


1.0 


2.43£18 


1£13 


15.2 


0.108 


2.23£27 


5.59£24 


399 


0.0044 


0.40 


1.0 


1£18 


1£13 


4.5 


0.334 


1.43E26 


2.13£23 


670 


0.052 


4.82 


1.0 


8£17 


1£13 


3.2 


0.423 


7.18£25 


9.66£22 


743 


0.095 


8.54 


1.0 


5£T7 


1E13 


1.4 


0.651 


1.65E25 


1.87E22 


885 


0.30 


26.93 


1.0 


2£17 


1£13 


0.1 


1.036 


6.92£23 


6.69£20 


1033 


2.21 


201.14 



TABLE IV: Q-ball properties and dark matter gravitino mass for n p = 10 and a„ = 0.01 as a function of varying M. 



m s { TeV) 


M( GeV) 


M,„(GeV) 


<P(0)/M,„ 


C0 2 /m 2 


E( GeV) 


Q 


E/Q( GeV) 


T d ( GeV) 


m 3/2 (GeV) 


1.0 


2.43£18 


5£13 


0.3 


0.940 


3.89£25 


5.88£22 


1004 


0.24 


20.14 


1.0 


2.43£18 


2£13 


4.4 


0.342 


5.33£26 


7.87£23 


667 


0.031 


2.55 


1.0 


2.43£18 


1.5£13 


7.4 


0.223 


9.23£26 


1.65E24 


558 


0.018 


1.30 


1.0 


2.43£18 


1£13 


15.2 


0.108 


2.23£27 


5.59£24 


399 


0.0044 


0.40 


1.0 


2.43£18 


9£12 


17.6 


0.092 


2.57£27 


6.93£24 


370 


0.0035 


0.32 


1.0 


2.43£18 


8£12 


25.8 


0.059 


5.31£27 


1.78E25 


298 


0.0016 


0.15 


1.0 


2.43£T8 


5£T2 


33.4 


0.042 


1.41E28 


4.01£25 


255 


0.0008 


0.076 


1.0 


2.43£18 


1£12 


502.4 


0.0008 


2.45£29 


6.71£27 


36 


3.49 x 10~ 6 


0.00031 



TABLE V: Q-ball properties and dark matter gravitino mass for n p = 10 and a a — 0.01 as a function of varying M, 



and 

nifermion = M m , (35) 

where A=< F$ > / < S > and M m = K < S >. This leads to gaugino masses for MSSM gauginos, A,,, at 1-loop 

mx i = n 5 A^- (36) 
47C 



and to scalar masses at 2-loops 



m? = 2n 5 A 2 £Q(£)- . (37) 

(=i 



- / a,- n 2 



471 



Here C, are the Casimir coefficients for the scalar multipet of interest, with i = 1,2,3 for the U (l)y, SU(2)l and SU (3) c gauge 
groups. The gravitino mass is related to < F$ > by the standard supergravity relation j26ll 

<Fs> „„ 

m3/2 = tit ' (38) 



therefore 



The scalar mass is then 



m i= i 



Therefore the theoretical gravitino mass is 



1 m s M m Mtc \ 

^v^^rfeJ' (41) 
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<p(0)/M,„ 


CQ 2 /m 2 


E( GeV) 


e 


E/Q 


Trf(GeV) 


m 3/2 (GeV) 


0.1 


15.22 


0.108 


2.23£27 


5.59£24 


398.87 


0.0044 


0.40 


0.01 


15.2 


0.108 


2.23£27 


5.59£24 


399.26 


0.0044 


0.40 


0.001 


15.19 


0.108 


2.23£27 


5.59£24 


399.51 


0.0045 


0.41 



TABLE VI: Q-ball properties and dark matter gravitino mass for n p = 10 as a function of varying a„. Here m s = 1 TeV, M,„ = 10 13 GeV and M = 2.43 x 
10 18 GeV. 



lip 


(p(0)/M m 


co 2 /m 2 


E{ GeV) 


Q 


E/Q( GeV) 


T rf (GeV) 


m 3/2( GeV) 


1 


38.0 


0.036 


2.23£28 


9.43£25 


236 


4.8£-4 


0.043 


5 


20.2 


0.079 


4.46£27 


1.30£25 


343 


2.4£ - 3 


0.21 


10 


15.2 


0.108 


2.23£27 


5.59£24 


399 


4.4£ - 3 


0.40 


25 


10.4 


0.162 


8.91£26 


1.85£24 


482 


LIE -2 


0.95 



TABLE VII: Q-ball properties and dark matter gravitino mass as a function of varying n p . Here m s = 1 TeV, M,„ = 10 13 GeV, M = 2.43 x 10 18 GeV and 
a o = 0.01. 



where VCa = (£? =1 Qa 2 ) 1 / 2 . 

For the case of the (u c d c d c ) 2 direction, VCoc is dominated by the SU(3) C contribution. C3 = 4/3 for an SU(3) C triplet and 
a3 m 1 /20 at the renormalization scale /j w Mm ~ 10 13 GeV. Therefore 

m 3/2 » 0.36 GeV x ( J-] ( ff» ) (J*-} f JLJ) . (42) 
3/2 \V"5*7 V 10 13 GeVy Vl TeV/ \20a 3 J 

Therefore with K ~ 1, values of m 3/ / 2 ~ 0.1 GeV can easily be achieved when M,„ ^ 10 13 GeV. We also note that K can be as 
large as 2-3 without violating perturbation theory, while «5 > 1 is also possible. 

This range of m 3 n and M m is consistent with the gravitino mass necessary for dark matter. For example, from Table 1 we 
find that with m s = 1 TeV, M — M p and M m = 10 13 GeV, the dark matter gravitino mass is W3/2 = 0.40 GeV. The theoretical 
gravitino mass in this case is W3/2 = 0.36 GeV/ ^/niK. So n$ = 1 and K ~ 1 would give a theoretical mass which is in agreement 
with the required gravitino dark matter mass. 



VII. CONCLUSIONS AND DISCUSSION 



We have shown that, under plausible assumptions regarding the parameters of the flat-direction potential and the evolution 
of the AD condensate fragments, both the baryon asymmetry and gravitino dark matter can originate from a d = 6 (u c d c d c ) 2 
flat direction of the MSSM with gauge-mediated SUSY breaking. The AD condensate fragments into Q-balls, which serve 
as a source of non-thermal stau NLSPs. These can subsequently decay to gravitinos without violating BBN constraints. The 
elliptical nature of the AD condensate in gauge-mediation plays a key role, leading to the formation of pairs of positive and 
negative global charge Q-balls, which enhances the gravitino number density and so lowers the gravitino mass below the 1 GeV 
upper limit from BBN. The number of gravitinos from Q-ball decay has an upper limit due to annihilation of the NLSPs. This 
prevents too many gravitinos being produced, which would otherwise heavily suppress the gravitino mass relative to the GeV 
scale, so reducing the messenger scale and resulting in Q-balls with too small E/Q to decay to NLSPs. 

There is a notable self-consistency to the model, d = 6 AD baryogenesis requires a low reheating temperature, Tr ~ 100 GeV, 
ruling out generation of gravitino dark matter via thermal scattering. Gravitino dark matter must then come from non-thermal 
NLSP decay (thermal relic NLSPs which decay to gravitinos being generally inconsistent with BBN). The Q-balls associated 
with d = 6 AD baryogenesis self-consistently provide the required source of non-thermal NLSPs as a result of their low decay 
temperature. Finally, the field strength at which d = 6 AD baryogenesis occurs is close to the messenger mass when OT3/2 is in 
the range 0.1-1 GeV, in which case Q-balls can form with E/Q large enough to allow them to decay to NLSPs and subsequently 
to gravitino dark matter while remaining consistent with BBN constraints. 

We have also shown that the gravitino mass which is necessary to account for the observed dark matter is easily consistent 
with the theoretical gravitino mass expected in a standard gauge-mediated SUSY breaking set-up. 

In the case of stau NLSPs, there is a preference for light staus, < 1 TeV, in order that is less than E /Q. This can be 
tested in the near future at the LHC 12411 . In addition, a relatively large messenger mass is predicted, in order that E/Q> itinlsp 
in the Q-balls. This large messenger mass scale may be detectable via the SUSY mass spectrum. In general, the stau and 
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sneutrino are favoured by BBN to be the NLSP in this model, since in that case m 3 / 2 can be close to 1 GeV. This should also be 
testable. Therefore, although it may be difficult to prove unambiguously that the model is correct, it would be significant if any 
of these signatures were observed experimentally. 

As well as direct experimental signatures, there may be more indirect signatures. The reheating temperature is predicted to be 
low, which may be observable via the gravity wave spectrum 12711 . In addition, AD baryogenesis can naturally produce baryon 
isocurvature perturbations, which may be observable [ 28, 29f] . 

Comparing with the thermal scattering model for gravitino dark matter production an obvious advantage of the AD 
baryogenesis/Q-ball decay model is that it simultaneously accounts for both gravitino dark matter and the baryon asymmetry. 
Since Affleck-Dine baryogenesis is a particularly simple mechanism for generating the baryon symmetry in SUSY models, the 
model presented here should be a strong candidate for the origin of gravitino dark matter in the gauge-mediated MSSM. As we 
have noted, the thermal scattering model and the AD baryogenesis/Q-ball decay model are mutually exclusive, corresponding to 
high (Tr ~ 10 7 GeV) and low (Tr ~ 100 GeV) reheating temperatures respectively. 

There are other models for gravitino dark matter based on Affleck-Dine baryogenesis. In |@] a model based on the early decay 
of relatively small charge Q-balls directly to non-thermal gravitinos was proposed. In 10] a model based on suppression of the 
Q-ball decay mode to NLSPs was proposed. The model presented here follows a more conventional approach, similar to the 
model of neutralino dark matter from Q-ball decay in the gravity-mediated MSSM I20I l2lll . 

The analysis we have presented here is very much a first study of the complete AD baryogenesis/Q-ball decay model for 
gravitino dark matter in the gauge-mediated MSSM. Our analysis makes a number of simplifying approximations. In particular, 
it uses a semi-analytical model for the initial fragmentation of the AD condensate and a simple phenomenological model for the 
evolution of the initial condensate fragments into an ensemble of positive and negative charged Q-balls. In reality, a distribution 
of Q-balls with different magnitude of charge would be expected, rather than the single charge assumed in the phenomenological 
model. In order to fully understand the evolution of the AD field and the final Q-ball distribution, a complete 3-D numerical 
simulation of the evolution of the GMSB AD condensate, similar to that conducted in II 1 Oil for the case of gravity-mediation, is 
necessary. The positive conclusion we have obtained here strongly motivates such a detailed numerical analysis. 

In conclusion, we have shown that that the baryon asymmetry and gravitino dark matter can originate from Q-ball decay in 
the gauge-mediated MSSM while remaining consistent with BBN constraints. In order to develop the model further, a detailed 
numerical analysis of Affleck-Dine condensate fragmentation and Q-ball decay is necessary. This will be the subject of future 
research. 
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Appendix: Homogenization and slowing of relativistic NLSPs from Q-ball Decay 

In this Appendix we discuss the slowing of the relativistic NLSPs from Q-ball decay via scattering from thermal background 
particles and the conditions for homogenization of the NLSP density. 

If we neglect the decay of the initial condensate fragments into pairs of ±Q-balls, at T c i the fragments would be separated by 
an average distance d s . In fact, the initial fragments are expected to decay to n p pairs of ±Q-balls, therefore the initial condensate 
fragments will be replaced by clusters of positive and negative Q-balls, in which the Q-balls are separated from each other by a 
distance d p at 7^. This is illustrated schematically in Figure 1. We will conservatively assume that d p <C d s , which will be true 
if the break up of the initial fragments into ±Q-ball pairs occurs well after the initial condensate fragmentation. (If the initial 
fragments immediately broke up into ±Q-ball pairs, the final separation of the Q-balls would be ~ d s /n p , therefore considering 
the separation to be d s is conservative.) In this case we need to understand whether the NLSPs from Q-ball decay can spread 
across the distance d s and become non-relativistic before annihilating. This is important because if they were unable to do so, 
the annihilation process would occur in inhomogeneous concentrations of NLSPs, in which the NLSPs would annihilate until 
the annihilation rate was less than the expansion rate. The resulting NLSP density would subsequently be diluted as the NLSPs 
spread out by thermal diffusion and homogenize. As a result, a much lower NLSP density would remain than would be expected 
under the conventional assumption of a homogenized non-relativistic NLSP density. This would then require a much larger 
gravitino mass to account for dark matter, violating BBN constraints on NLSP decay. 

We first consider the evolution of NLSPs from Q-ball decay in the absence of annihilations. The distance d s can be estimated 
from the initial fragmentation of the AD condensate. The initial condensate fragments have an initial diameter X f rag and therefore 
their centres are initially separated from each other by a distance approximately A. /«,<,. These fragments then expand from each 
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other as the Universe expands. Their separation when the Q-balls decay at T d is therefore 

!/ 3 T_ / U, \ 2 / 3 



ds fl(7fi) a fmg " Jrag \g{T d ) ) T d \H(T R )J 

Once a relativistic NLSP is produced from the decay of a Q-ball, there will be two stages to its subsequent motion. Firstly, 
it will stream out at close to the speed of light until scattering from thermal background particles slows it to a non-relativistic 
velocity, which we define to occur after a distance cInr- The NLSP will then thermalize and diffuse out from the region where 
the NLSPs become non-relativistic. 

If cLnr > ds then the NLSPs from decay of the ensemble of Q-balls will homogenize by relativistic free-streaming before 
becoming non-relativistic. On the other hand, if d^R < ds, we would have to consider the time for the non-relativistic NLSPs to 
thermally diffuse out over the distance between the Q-balls. In practice, we find that d^R 3> ds in the cases of interest to us here, 
therefore diffusion will play no essential role. 

We define At^R to be the time for the NLSP to become non-relativistic. We will focus on the case of stau NLSPs. Since staus 
are electrically charged, they will scatter at tree-level from thermal background photons. (The scattering rate of sneutrino and 
other neutral NLSPs from thermal background particles will be generally smaller, in which case AInr will be longer. Therefore 
dpjR for staus can be considered a lower bound on d^R for NLSPs in general.) The scattering cross-section for relativistic staus 
from thermal photons is 

(A-2) 
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This assumes that the energy E of the staus satisfies \2ET < m\, where the mean energy of the thermal photons is ps 3T. The 
scattering rate of the relativistic staus is then F sc = riyOv, where riy = 2AT 3 /n 2 is the background photon density and the relative 
velocity is v — 1 . Therefore 

1.2 e 4 r 3 
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The average energy loss per scattering is 



AE = 2 [ x 3T . (A-4) 
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Therefore the energy decreases to E fx in time 

(x-l)ml , 

At = - -—TZ 1 . (A-5) 

6 ET sc 

We will define the time to become non-relativistic At^R to correspond to E = Ei/2, where £, is the initial stau energy, in which 
case x — 2 and 



1 mi 



2 



Ar-=6#r- (A-6) 

(To be more precise we should consider x = Ei/ m x , but this will give a value for AtpjR which differs from Eq. ( IA-6I > by only a 
small factor.) The distance travelled by the stau NLSPs before becoming non-relativistic is therefore 

1 4 
7t lfl~ 

d NR fa cAt NR = — ^r^4 ■ (A-7) 

If dNR > ds then the staus from Q-ball decay will become homogeneous via relativistic free-streaming. In this case there is 
no need to consider the subsequent thermal diffusion of the non-relativistic NLSPs. From Table 3 we find that this condition is 
easily satisfied in practice. 

We also need to check that the injected NLSPs become non-relativistic before annihilating. We therefore compare the 
mescale At^R with the timescale for relativistic stai 
tivistic staus (via t-channel x exchange to photons) is 



(A-8) 
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where E 2 ^> to?. The number density of staus will generally be extremely small compared with the number density of background 
photons, since even if we consider all the Q-balls to instantaneously decay to staus at T c / without annihilation (which greatly 
exaggerates the number density of relativistic NLSPs), the total NLSP density would satisfy n^LSp/ny ~ O(l00)riB/riy ~ 10 -8 , 
where we have assumed that the energy in the initial fragments is at most 0(100) times the baryonic charge times the AD scalar 
mass. The annihilation cross-section Eq. dA-8b is also less than the scattering cross-section with background photons Eq. ( 1A-21 >. 
The maximum annihilation rate for relativistic staus from Q-ball decay therefore satisfies 

Vamrel = iNLSPGannrel < ^^~ % nfS amre l < lO~ & tl y O sc = 10~ 8 r je . (A-9) 

Therefore At ann re \ ^> lO 8 ^. 1 . Since, from Eq. (IA-6I 1. At^R is at most ~ lO 4 ^. 1 (using m% ~ 500 GeV and T ( i ~ 5 MeV), we find 
that At ann K i ^> 10 4 A?tvr. So, rather generally, the relativistic NLSPs will become non-relativistic long before they can annihilate. 
In this case there will be annihilation of a homogeneous density of non-relativistic stau NLSPs, as in the conventional Q-ball 
decay scenario. 

Thus if cLnr > ds, the conventional picture of a homogeneous non-relativistic NLSP density annihilating to a final non-thermal 
NLSP relic density will apply. 
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